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Abstract 

We propose a parametrization of two-dimensional geophysical turbulence in the form of a re- 
laxation equation similar to a generalized Fokker-Planck equation (Chavanis 2003). This equation 
conserves circulation and energy and increases a generalized entropy functional determined by a 
prior vorticity distribution fixed by small-scale forcing (Ellis et al. 2002). We discuss applications 
of this formalism to jovian atmosphere and Jupiter's great red spot. We show that, in the limit of 
small Rossby radius where the interaction becomes short-range, our relaxation equation becomes 
similar to the Cahn-Hilliard equation describing phase ordering kinetics. This strengthens the 
analogy between the jet structure of the great red spot and a "domain wall". Our relaxation equa- 
tion can also serve as a numerical algorithm to construct arbitrary nonlinearly dynamically stable 
stationary solutions of the 2D Euler equation. These solutions can represent jets and vortices that 
emerge in 2D turbulent flows as a result of violent relaxation. Due to incomplete relaxation, the 
statistical prediction may fail and the system can settle on a stationary solution of the 2D Euler 
equation which is not the most mixed state. In that case, it can be useful to construct more gen- 
eral nonlinearly dynamically stable stationary solutions of the 2D Euler equation in an attempt to 
reproduce observed phenomena. 
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I. INTRODUCTION 



An important problem in geophysical fluid turbulence is to obtain a small-scale 
parametrization of the 2D Euler equation in order to compute only large-scale motion. In- 
deed, in the limit of vanishing molecular viscosity (which is a limit relevant for geophysical 
flows), the 2D Euler equation develops a complicated filamentation at smaller and smaller 
scales. This is a redoutable problem for numerical simulations since those filaments rapidly 
reach the mesh size and generate numerical instabilities. Therefore, one usually adds an 
"artificial" or "numerical" viscosity > in the Euler equation in order to smooth-out the 
small scales and regularize the problem. This viscosity is a simple parametrization of the 
turbulent Reynolds tensor which accounts for the correlations of the velocity fluctuations. 
However, this artificial viscosity breaks the conservation laws of the inviscid dynamics, in 
particular the conservation of energy, and leads ultimately to a state of rest while it is known 
that 2D flows generate long-lived coherent structures. One can in principle limit these prob- 
lems by using a small viscosity but this is computationaly costly because we have to solve 
the small-scales in order to solve the large scales which are of physical interest. 

How can we improve the Reynolds parametrization and what should we add to diffusion 
in order to restore the conservation of energy and obtain stable stationary states on the 
coarse-grained scale? It is illuminating to realize that this problematic is related to the one 
encountered by Einstein in his investigation of Brownian motion and by Chandrasekhar in 
his stochastic approach of stellar dynamics [l( . The solution is now well-known: they solved 
the problem by introducing a dynamical friction in addition to diffusion so as to recover the 
correct equilibrium state. In their cases, the equilibrium state is the Maxwell-Boltzmann 
distribution predicted by statistical mechanics. The condition that the Maxwell-Boltzmann 
distribution is the stationary solution of the Fokker-Planck equation leads to the celebrated 
Einstein relation £ = Dj3 between the coefficient of diffusion and the friction. Can we follow 
a similar approach in 2D turbulence and which adaptations should we make? 

Two theories have been developed to predict the (meta)-equilibrium state resulting from 
violent relaxation in terms of statistical mechanics. The approach of Miller Q] and Robert 
& Sommeria 3| applies to unforced systems which are strictly described by the 2D Euler 
equation. It predicts the most probable distribution of the flow depending on the energy, the 
circulation and on the infinite set of Casimirs (or vorticity moments) that are conserved by 
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the Euler equation. This approach is quite powerful in the case of simple initial conditions 
made of a small number of patches for which the number of control parameters is reduced. It 
may also describe more general situations in the limit of strong mixing or low energy studied 
by Chavanis & Sommeria j^J which is characterized by a finite number of relevant contraints 
(the first moments of vorticity). However, the approach of also presents limita- 

tions that have been discussed by several authors. In particular, the statistical prediction 
presupposes that the (fine-grained) initial conditions are known. Indeed, the equilibrium 
state depends not only on circulation (the equivalent of mass) and energy as in ordinary 
statistical mechanics but also on the values of all the vorticity moments T n = j u n d 2 r which 
cannot be determined from the coarse-grained field when the vorticity has mixed (for t > 0) 
since uj n 7^ uf 1 . Since the initial conditions are not known in general (and we are never sure 
whether a flow is already mixed or not), we cannot really make statistical predictions. This 
approach also assumes that the evolution is inviscid and freely evolving. However, in real 
situations this is not the case and the conservation of all the Casimirs is abusive. 

An alternative approach has been developped by Turkington and collaborators j^, |f| in 
order to remedy these difficulties. This approach is particularly relevant for geophysical 
and astrophysical flows that are forced at small-scales. For example, convective plumes are 
known to form regularly in the atmosphere of Jupiter. Their significant horizontal divergence 
coupled with rotation leads to small-scale anticyclonic vorticity production. Therefore, there 
is a permanent forcing and dissipation at small scales. These effects destroy the conservation 
of the Casimirs while establishing instead a global distribution of vorticity. In that case, it 
may be more relevant to fix a prior distribution of vorticity x(°~) instead of the Casimirs. This 
amounts to treating the constraints associated with the fragile integrals T n > 2 canonically 
while the robust integrals V and E are treated microcanonically. In this point of view, we 
abandon the strict conservation of all the Casimirs and keep their influence indirectly in 
the function y(cx). This prior distribution has to be adapted to the geophysical situation as 
discussed in [6]. This is similar to our idea [7] of fixing a generalized entropy S = — f C(q)d 2 r 
a priori, depending on the situation contemplated. In fact, we will show that the prior 
vorticity distribution x(°") directly determines a form of generalized entropy C(q). 

Once a prior vorticity distribution is given, the generalized entropy is known and the 
equilibrium state depends only on the value of circulation and energy (robust constraints). 
We are thus in a situation of usual thermodynamics but with a generalized form of entropy. It 
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is therefore possible to apply phenomenological arguments similar to those given by Einstein 
and Chandrasekhar to derive a relaxation equation towards the equilibrium state which 
conserves energy and circulation and increases the generalized entropy. This equation has the 
form of a generalized Fokker-Planck equation [7j involving a turbulent viscosity (diffusion) 
and a drift. As in Brownian theory, the diffusion and the drift coefficients are connected by 
a sort of Einstein relation involving here a negative temperature. 

In this paper, after a brief summary of the two statistical theories of 2D turbulence 
(Sec. [n}, we give in Sec. IIHI a more precise construction and physical justification of the 
relaxation equation introduced in [7], emphasizing the connexion with the work of Ellis et 
al. 0. We present a new derivation of this generalized Fokker-Planck equation by using 
an analogy with Onsager's linear thermodynamics and with the theory of Brownian motion. 
The domain of application of our parametrization is therefore better established and the 
meaning of the generalized entropy, determined by the prior vorticity distribution, clearer. 
For given prior vorticity distribution xi a )i there is no free parameter in our model. In Sec. 
IIV| we apply this formalism to the atmosphere of Jupiter. We consider the limit of small 
Rossby deformation radius, corresponding to short-range interactions, and show that our 
relaxation equation reduces to a form of Cahn-Hilliard equation that is used in the theory 
of phase ordering kinetics [8]. Its stationary solutions describe "domain walls" accounting 
for the annular jet structure of Jupiter's great red spot. On the other hand, according to 
the refined stability theorem of Ellis et al. [6], our relaxation equation can also be used 
as a numerical algorithm to construct arbitrary nonlinearly dynamically stable stationary 
solutions of the 2D Euler equation. This can be used to construct confined structures 
that are not described by the statistical theory because of incomplete relaxation (lack of 
mixing/ergodicity) [j]]. This paper discusses the main ideas and presents the formalism. 
Numerical simulations illustrating these results will be given in a future work. 

II. SUMMARY OF THE STATISTICAL THEORIES 
A. The quasi-geostrophic equations 

Large-scale oceanic and atmospheric flows are quasi two-dimensional due to planetary 
rotation and stratification 10]. This property prevents vortex stretching and can lead to 



4 



the formation of large-scale coherent structures (e.g., Jupiter's great red spot, cyclonic and 
anticyclonic atmospheric eddies,...) with robust stability (Jupiter's great red spot was ob- 
served more than three centuries ago). In addition, in geophysical situations, the Reynolds 
numbers are so high that viscosity is not expected to play a crucial role for the timescales 
and length scales of interest. A relevant model is provided by the quasi-geostrophic (Q.G.) 
equations 

| + u-Vg = 0, (1) 

q = -Ai/j + ^-h(y), (2) 

where q is the potential vorticity (PV), i/j the stream-function defined by u = — z x Vip, R is 
the Rossby radius of deformation and h(y) is an effective topography including the /3-effect. 
The second term in Eq. (J2J) is similar to the Debye shielding in plasma physics. More 
generally, we can consider a relation of the form ip(r,t) = j u(\r — r'\)[q(r',t) + h(y')]d 2 r' 
where u(|r — r'|) is a binary potential of interaction. In the Q.G. model, one has w(£) = 
(1/27t)K q (^/R), where K Q (z) is a modified Bessel function. It has the asymptotic behaviours: 
K (z) ~ — In z for z — > and Kq(z) ~ \p^ e ~ z f° r z ~~ * +oo. Stationary solutions of the 
Q.G. equations are specified by a relation of the type q = f(ip) where / is an arbitrary 
function. 

The Q.G. equation (JTJ) simply states that, in the absence of viscosity, the potential vor- 
ticity q is conserved by the flow. The Q.G. equations conserve the energy 

E = \ f(q + h)1><pT, (3) 

and the Casimirs Ch = J h(q)d 2 r, where h is any continuous function of q. This infinite 
set of constraints is due to the inviscid dynamics and results from the conservation of q 
and the incompressibility of the flow. The conservation of the Casimirs is equivalent to the 
conservation of the PV moments 




which include the circulation T and the enstrophy T2- The moments T n with n > 2 will 
be called fragile constraints because they are affected by viscosity or coarse-graining (since 
q n 7^ q n ) while the circulation V and energy E will be called robust constraints because 
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they are approximately conserved on the coarse-grained scale or in the presence of a small 
viscosity. 



B. The statistical equilibrium state 

The Q.G. equations generate a complicated mixing process and develop filaments at 
smaller and smaller scales. If we work with a finite resolution e, the coarse-grained PV 
q(r,t) (the local average of q(r,t) over a cell of size e 2 ) will reach a metaequilibrium state 
on a very short timescale. This is an inviscid violent relaxation driven by purely inertial 
effects (mixing). In practice, a natural coarse-graining is played by inherent viscosity which 
smoothes out the small scales. The resulting structure has the form of a large-scale vortex 
or a jet. 

There has been some attempts to describe this metaequilibrium state in terms of sta- 
istical mechanics P, 0, Q| . The same ideas had been developed previously by Lynden-Bell 
ri| in stellar dynamics for the Vlasov-Poisson system (see [3, 3 for a description of this 
analogy). Assuming that mixing is efficient, which is linked to a condition of ergodicity, the 
most probable local distribution of vorticity p(r, a) is obtained by maximizing the mixing 
entropy 

S[p] — — I plnp d 2 rda, (5) 



taking into account all the constraints imposed by the dynamics. This leads to the Gibbs 
state 

P(r ) a) = ^x(^)e- (W+Q)CT , (6) 

where x( a ) = ex P(~ E n >i ^n 17 ") accounts for the conservation of the fragile moments 
r ra >i = / pcr n d 2 rdcr and the "partition function" Z = J x(cr)e~^^ +a - )cr (io" is determined by 
the local normalization condition J pdo = 1 (a and [3 are the usual Lagrange multipliers as- 
sociated with the robust constraints T and E). Using the Gibbs state (jBJ), the coarse-grained 
PV q = J pada is related to the streamfunction according to 

~ q = -)rw = F w+ a ) = fw- ( ? ) 

Thus, for given initial conditions, the statistical theory selects a particular stationary solution 
of the Q.G. equations (on the coarse-grained scale) among all possible ones. Specifically, the 
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equilibrium state is obtained by solving the differential equation 

-&l> + -^l>-Kv) = f anJl (il>), (8) 

and relating the Lagrange multipliers a n , f3 to the constraints T n , E. In this approach, the 
function x{°~) is determined from the constraints a posteriori. Indeed, we need to solve the 
full problem in order to get the expression of x{°)- in this sense, the vorticity moments are 
treated microcanonically. Taking the derivative of Eq. (|7|) ; it is easy to show that 

q'W = -PQ2, q 2 = J p(a - qfda > 0, (9) 

where q<i is the local centered enstrophy. Therefore, f(ip) is a monotonic function increasing 
at negative temperatures and decreasing at positive temperatures. 



C. Prior vorticity distribution and generalized entropy 

For freely evolving 2D turbulence, the function x{ a ) is determined by the initial conditions 
through the values of the Casimir integrals. However, in geophysics, there exists situations in 
which the flow is continuously forced at small-scales so that the conservation of the Casimirs 
is destroyed. Ellis et al. |6j have proposed to take into account these situations by treating 
the function x( a ) canonically (see Appendix EJ. Therefore, the function x( a ) must be 
viewed as a prior vorticity distribution fixed by the small-scale forcing. Its specific form 
has to be adapted to the situation. In that context, the local vorticity distribution can be 
written p(r, a) = x(cr)yu(r, a) where xi a ) is given. The optimal distribution is then obtained 
by maximizing the relative entropy [6j]: 

SM=-/Mr,»)h,Mr,„)x(«r)rfW, (10) 

at fixed mass, energy and normalization (no other constraint). This leads to Eq. © but with 
a different interpretation. The local distribution of vorticity p(r, a) is equal to the global 
vorticity distribution x( a ) locally modified by the large-scale flow through the Boltzmann 
factor ^-^e~^^ +a ^ a . The coarse-grained vorticity is still given by Eq. (jZJ. However, the 
crucial point to realize is that the function F is now entirely determined by the prior vorticity 
distribution x( a ) while, in the former case, it had to be related to the initial conditions in a 
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very complicated way. This makes this new approach very attractive. The resulting macro- 
distribution q(e) = J x(o")e _cre (io" with e = + a is formally similar to the notion of 
superstatistics developed by Beck & Cohen [ljj], see [3] for the description of this analogy. 

We now show that the prior distribution function determines a form of generalized en- 
tropy. Since the relation q = F(f3ip + a) is monotonic, it can be obtained by maximizing a 
functional of the form 

S = - [ C(q)d 2 r, (11) 



where C is a convex function (i.e. C" > 0), at fixed circulation and energy. Indeed, 
introducing appropriate Lagrange multipliers and writing the variational principle in the 
form 

5S - (35E - a5T = 0, (12) 
we find that the critical points of S at fixed circulation and energy are given by 

C\q) = -j3i) - a. (13) 

Since C is a monotonically increasing function of q, we can inverse this relation to obtain 

q = F(M + a) = f(i>), (14) 

where 

F(x) = (C'y'i-x). (15) 

From the identity 

= -P/C"(q), (16) 

resulting from Eq. (jlHj) . q(i/j) is monotonically decreasing if f3 > and monotonically increas- 
ing if (3 < 0. Thus Eq. is equivalent to Eq. (jZJ provided that we make the identification 
()15|). Therefore, for any function F(x) determined by the prior vorticity distribution x(cr), 
one can associate a generalized entropy (fTTj) where C(q) is given by Eq. (fTB]) or equivalently 

C(5 ) = -/V(^. (17) 

Furthermore, S\g] really has the status of an entropy in the sense of large deviations. Indeed, 
Ellis et al. jfij show that the distribution probability of the coarse-grained vorticity field can 
be written in the form of Cramer's formula P\q] ~ e n5 ® where n is the number of sites of 



S 



the underlying lattice introduced in their mathematical analysis. Therefore, the optimal PV 
field q maximizes S[q] at fixed circulation and energy. They also relate the entropy S[q] to 
the prior distribution x( a ) by using the formalism of large deviations. The steps pijl -(|17 )l 
presented above are a simple alternative to construct the generalized entropy associated with 
It of course leads to the same result as p. 

D. Refined condition of nonlinear dynamical stability 

Ellis et al. also show that the microcanonical variational principle 

Max S[q] | E[q] = E, T[q] = T, (18) 

provides a refined condition of nonlinear dynamical stability for the Q.G. equations. It is 
stronger than the canonical variational principle 

Max J[q] = S[q]-/3E[q] | T[q\ = T, (19) 

which just provides a sufficient condition of nonlinear dynamical stability (J is equal to the 
free energy F = E — TS multiplied by — (3 where /3 = 1/T is usually negative). The two 
variational principles differ if the ensembles are not equivalent which is the case for systems 
with long-range interactions. This can have important implications for geophysical flows 
since vortex solutions can be stable in the microcanonical ensemble while they are unstable 



in the canonical ensemble 



a 



This is similar to "core-halo" structures with negative 
specific heats in astrophysics |17|. Note that, in the view point of Ellis et al, the statistical 
equilibrium macrostates are stable with respect to the perturbations on the microscopic 
scale (thermodynamical stability) and the steady mean state q is stable with respect to 
macroscopic perturbations (nonlinear dynamical stability). 

It is worth mentioning, however, that the results of Ellis et al. [6] on nonlinear dynamical 
stability are valid independently of the statistical theory. They apply to any stationary 
solution of the Q.G. equation, whether or not it corresponds to the statistical equilibrium 
state (J7J). This is important in the case of incomplete relaxation where the statistica 
diction fails 9]. The stability results still apply but in that case S is a H-function 
not an entropy and J is an Energy-Casimir functional [3], not a free energy. Furthermore, 
the variational principles ()18j) and (|19j) have not a direct thermodynamical interpretation 
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in that case. However, because of the formal ressemblance of these variational principles 
with those of thermodynamics, it is possible to develop a thermodynamical analogy [7( to 
settle the nonlinear dynamical stability of a stationary solution of the Q.G. equations. In 
this context, the "canonical" stability criterion (J19|) is equivalent to the criterion of formal 
nonlinear stability given by Holm et al. |l8j . The Arnold's theorems provide particular 
conditions of formal nonlinear stability in the sense of ()19|). The "microcanonical" stability 
criterion (fTHj) proved by Ellis et al. jy] refines these theorems. In particular, a flow can be 
stable by the criterion ()18|) although it violates the Arnold's theorems. 

We note finally that the stability criterion (jl8)l is consistent with the phenomenology of 2D 
turbulence. Indeed, the H-functions increase on the coarse-grained scale (in the sense that 
S'IMO] — ^[^(O)] where g(0) = g(0) is un-mixed Q|) or in the presence of a small viscosity 
(since S = v j C"(q)(Vq) 2 d 2 r > 0), while the energy and the circulation are approximately 
conserved. Therefore, the metaequilibrium state is expected to maximize a certain H- 
function at fixed circulation and energy. This generalizes the concept of "selective decay" 
which considers the decrease of enstrophy at fixed energy and circulation. However, the 2D 
Euler equation does not select a universal if -function. Minus the enstrophy — T2 = — j q 2 d 2 r 
and Tsallis Q-functionals Sq = — ^pi j (,Q®~ q)d 2 r are particular if-functions (not entropies) 
that sometimes occur in 2D turbulence in case of incomplete relaxation as in the plasma 

n 

experiment of Huang & Driscoll [19|- They are associated with confined stationary solutions 
of the 2D Euler equation that have a compact support ("polytropic vortices") 0. Note 
that many other confined solutions can emerge in 2D turbulence, e.g. [^j, which are not 
described by Tsallis of by the enstrophy functionals. 

III. A PARAMETRIZATION OF GEOPHYSICAL FLOWS 

In the context of freely evolving 2D turbulence, a thermodynamical parametrization of 
the 2D Euler equation has been proposed by Robert & Sommeria (^J in terms of relaxation 
equations based on a Maximum Entropy Production Principle. These equations conserve 
all the Casimirs, increase the mixing entropy @ and relax towards the Gibbs state ©. 
In the situation considered by Ellis et al. p] where the system is forced at small-scales, 
we can propose an alternative parametrization of the 2D Euler equation. In that case, we 
have seen that only the energy and the circulations are conserved. The conservation of the 
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Casimirs is replaced by a prior vorticity distribution x( a )- This fixes the form of generalized 
entropy. We can thus propose to describe the large-scale evolution of the flow in that case by 
a relaxation equation which conserves energy and circulation and increases the generalized 
entropy (jTTjl associated with x(er). 

Taking the local average of the Q.G. equation JI}, we have 

§ + u-Vg = -V-J, (20) 

where J = qu is the Reynold's diffusion current. We shall determine its expression by 
using a procedure similar to Onsager's linear thermodynamics. Noting that the "chemical 
potential" a(r,t) = —f3ip — C'(q) is uniform at equilibrium according to Eq. (fT3j). we take 
J = D'Va. This choice yields a parametrization of the form 



(21) 



where, for convenience, we have droped the bars and set D = D'C"(q). The convex func- 
tion C (q) is a non-universal function which encapsulates the complexity of the fine-grained 
dynamics and depends on the situation contemplated. It is determined by the prior distribu- 
tion x( cr ) which encodes the statistics of vorticity fluctuations induced by small-scale forcing 
Equation (J2T|) can be interpreted as a generalized Fokker-Planck equation j^. We note 
that the term arising in the eddy flux in addition to the usual diffusion is a systematic drift 
— £V?/>. This drift must exist in order to yield the correct distribution at equilibrium. The 
necessity of this drift is similar to the necessity of the dynamical friction in stellar dynamics 
(see the first, phenomenological, approach of Chandrasekhar |l] on dynamical friction). The 
drift coefficient (mobility) is given by £ = D/3 which we link to the celebrated Einstein re- 
lation. Thus, the systematic drift is the counterpart of the ordinary dynamical friction and 
the preceding relation is the counterpart of the "fluctuation-dissipation" theorem in Brow- 
nian theory. This drift is usually absent from usual modelling of 2D turbulence where only 
a viscosity (or hyperviscosity) is introduced. This corresponds to the infinite temperature 



limit {(3 = 0) of our parametrization. 
from kinetic theories as attempted in 



lie physical origin of the drift should be understood 
] using a quasilinear theory of the Euler equation 
or in 2'_| for point vortices. Unfortunately, these formal methods do not describe the very 
nonlinear phase of "violent relaxation" so that more phenomenological methods are prefered 
in that case. 
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In Onsager's linear thermodynamics, the inverse temperature (3 is constant. This is 
appropriate to a canonical situation. In that case, the generalized Fokker-Planck equation 
()21|) increases the free energy J = S — (3E at fixed circulation provided that the diffusion 
coefficient D is positive. However, statistical ensembles can be inequivalent in 2D turbulence 
and the correct description is the microcanonical ensemble. The conservation of energy 
(J3J) can be taken into account by determining the evolution of f3(t) in Eq. (}2"Tj) from the 
requirement E = as in 2^. This yields 

m = J DV %I^\ (22) 

With this modification, the generalized Fokker-Planck equation (J2"T|) increases the gen- 
eralized entropy S at fixed circulation and energy until the maximum entropy state is 
reached. It automatically determines the right values of the Lagrange multipliers a and 
(3 = lim t ^ +00 (3{t) at equilibrium. This relaxation equation selects only maxima of entropy, 
not minima or saddle points that are linearly unstable Q. Note that the relaxation equa- 
tion (12 lj) with (J22|) can also be obtained from a Maximum Entropy Production Principle, 
see |7|, by using a procedure similar to that employed by Robert & Sommeria 23] in their 
parametrization. 

Heuristic arguments Q| or more formal kinetic theory suggest that the turbulent 
diffusion in Eq. (J21|) is given by D = Ke 2 q^ 2 where K is a constant of order unity, e is the 
scale of unresolved fluctuations and q 2 the local centered enstrophy (see Appendix \B§ . This 
space dependant diffusion coefficient, related to the fluctuations of vorticity, can account for 
incomplete relaxation 2M an d was used in oceanic modelling ^25] . Now, from Eqs. (|T^|l and 



we have the important relation 



* = c4r (23) 



which can be seen as a sort of fluctuation-dissipation theorem |26[ . This equation is valid at 
equilibrium but, in the spirit of linear thermodynamics, we shall still use it out of equilibrium 
as a convenient approximation (see also Appendix [HJ) • In that case, the diffusion coefficient 
can be expressed in terms of q alone as 

D = ^. (24) 
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Note that we can also obtain the parametrization (|2ip by closing the hierarchy of moment 
equations of j^jj with relation (j2HJ), see Q. In this sense, (J25|) can be regarded as a closure 
relationship. 

We have thus obtained a parametrization of 2D geophysical flows with no free parameter 
once the prior distribution of vorticity is fixed. These equations are expected to be valid 
close to the equilibrium state in the spirit of Onsager's linear thermodynamics. However, 
they may offer a useful parametrization even if we are far from equilibrium. Alternatively, 
according to the refined stability criterion of Ellis et al. (| , the relaxation equation ()21|) can 
also provide a poweful numerical algorithm to compute arbitrary nonlinearly dynamically 
stable stationary solutions of the Q.G. equations characterized by a function q = f(ip) where 
/ is related to C by Eq. ()15|). This can be used to construct confined solutions that are 
not consistent with statistical theory because of incomplete relaxation. We need just start 
from an initial condition (guess) with given (E, T) and Eqs. ()21j ) -([22| ) will relax towards a 
stationary solution specified by / or C which is, by construction, a nonlinearly dynamically 
stable stationary solution of the Q.G. equations with these prescribed values of energy and 
circulation. In that context, since we are only interested by the equilibrium state, we can 
take D = Cst. and even drop the advective term in Eq. (j21j) . Another numerical algorithm 
has been proposed by Turkington & Whitaker |27( based on iterative schemes techniques and 
it was used with success in many situations. The relaxation equations (|21 )) -([22 p can offer an 
alternative numerical algorithm. We again emphasize that constructing exact, nonlinearly 
dynamically stable stationary solutions of the 2D Euler equation is an important problem 
in itself and is an interest of the relaxation equation (J2T|) . independently of the statistical 
theory. In particular, it is difficult to solve the differential equation (JHJ) directly (for flows 
without simple symmetries) and be sure that the solution is stable. 

Some examples of prior vorticity distributions and corresponding generalized entropies 
and parametrizations have been collected in For example, in the two- levels case, the 

generalized entropy is similar to the Fermi-Dirac entropy 



where q = po\ + (1 — p)ao- We note that, in this particular case, the generalized entropy 
S[q) coincides with the mixing entropy S[p] defined in (J5J) since p(r, a) = p(r)S(a — a{) + 
(1 — p(r))8(a — do) where p is expressed in terms of q. It also coincides with the relative 




(25) 
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entropy S[p] denned in (|10|) since %(cr) = 5(cr — a\) + — °o)- This implies that the 
two-levels case can be seen either as the statistical equilibrium state resulting from the free 
merging of large patches with two types of vorticity (point of view of Sec. I11BJ) or as a 
particular type of prior vorticity distribution corresponding to a small-scale forcing with 
two intense peaks of vorticity (point of view of Sec. Ill C|) . Although the equations are the 
same, their interpretation is different. We think that the second interpretation is the most 
relevant for geophysical flows due to the existence of convective plumes. On the other hand, 
in the geophysical situation considered by Ellis et al. 6], the prior vorticity is taken 
decentered gamma function and the generalized entropy is 



cLS cL 



q - ~ ln(l + Xq) 
A 



d 2 r, (26) 



1 /2 

where Q 2 and 2\Q 2 are the variance and the kurtosis of x( a )- F° r A — > 0, S[q] is propor- 
tional to minus the enstrophy T 2 = j q 2 d 2 r. Note that minus the enstrophy can represent 
either a H-runction arising ie case of incomplete relaxation as in Q, a generalized entropy 
associated with a gaussian distribution x( <T ) as m |2j or the approximate form of mixing 
entropy © in the limit of strong mixing as in |4| . It can thus have several interpretations. 

These results can be readily extended to the more realistic shallow-water (SW) equations 
whose statistical mechanics has been developed by Chavanis & Sommeria |26j in the case of 
unforced flows. To apply the point of view of Ellis et al. we must write the total vorticity 
distribution as p(r, cr) = x(c")M r , <T ) where x( CT ) is given. The optimal distribution p*(r, a) 
is then obtained by maximizing the entropy S[p] = — j fi(r, a) ln/i(r, a)x{o~)h{r)d 2 Ydij (h 
is the elevation of the fluid) at fixed mass, circulation and energy (no other constraints). 
This leads to the same results as in 2^j], equivalent to Eq. ©, except that now the prior 
vorticity distribution x(°") must be regarded as given. The coarse-grained PV maximizes a 
generalized entropy S[q] = — J C(q)hd 2 r at fixed energy, mass and circulation, where C(q) 
is determined by x( a ) as i n Sec. Ill CI Relaxation equations equivalent to Eq. (|2*Tj) are 
obtained by maximizing S[q] at fixed E. This yields the same relaxation equations as in [26j 
except that now q 2 = 1/C"(q). 
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IV. THE LIMIT OF SMALL ROSSBY RADIUS 



We shall now discuss some applications of these ideas to geophysical and astrophysical 
flows. The application of the statistical theory of Miller-Robert-Sommeria to jovian 

flows and Jupiter's great red spot has been developed in j3, 29, 30 , [jlj]. In particular, 
Bouchet & Sommeria |3(| consider a limit of small deformation radius R — > and interprete 
the annular jet structure of Jupiter's great red spot (GRS) as the coexistence of two thermo- 
dynamical phases in contact, like in the van der Waals phase transition (an analogy sketched 
in 2sl). They show in particular that a very good agreement with the GRS can be obtained 
in the two-levels approximation of the statistical theory, and that this approximation can be 
motivated by the existence of convecting plumes in the jovian atmosphere. Turkington and 
collaborators ^, |(| have also developed a model of Jupiter's great red spot based on their 
statistical mechanics approach. In particular, they consider a realistic topography (Limaye 
profile) and show that a vortex structure appears precisely at the latitude of Jupiter's great 
red spot. They also emphasize the inequivalence of statistical ensembles and the relevance 
of a microcanonical formulation to settle the stability of jets and vortices. 

We shall here discuss some connexion between the relaxation equation (|21j) and the Cahn- 
Hilliard equation of phase ordering kinetics. This connexion will strengthen the resemblance 
of Jupiter's great red spot with a "domain wall" . The relaxation equation (J2~T|) can be written 



dq 
dt 



+ u • Vg = -V • 



D 6J 



C»{q) 5q 



(27) 



where 5/Sq is the functional derivative and where we have introduced the free energy 

J[q] = - J C(q)d 2 r - ~/3 J ' (q + h)^d 2 r. (28) 

We shall now consider the limit of small deformation radius R — > leading to short-range 
interactions. This limit is appropriate to jovian vortices such as the great red spot (GRS). 
In the limit of small deformation radius — ^ 0, the relation (j2J) between q and ip can be 
expanded to second order in R 2 as 



ip ~ R 2 (q + h) + R 4 Aq. 



(29) 



To simplify the discussion further, we shall assume that C(q) is an even function of <fi = q—B, 
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where B is a constant. Then, the free energy (|28|) can be written 

J[<f>] = (3R 2 j \^R 2 (S7<P) 2 + V(<f>) - 0/ijd 2 r, (30) 

with V(</>) = —C(4>) I (5R 2 — (l/2)0 2 (we have used the freedom to introduce a term Xq in 
the free energy to eliminate the constant B). In this limit of short-range interactions, the 
free energy is equivalent to the Landau functional (times PR 2 ). Since (3 is negative in cases 
of physical interest, we have to minimize the Landau functional. Note that for short range 
interactions, the ensembles are equivalent so that we can directly work with the free energy 
functional (|3U}l instead of the entropy. In addition, the relaxation equation (|2*T|) becomes 
similar to the Cahn-Hilliard equation 

f = pR2v ' { cW) V(i?2A0 " V ' {4>) + h) }' (31) 



which appears in the theory of phase ordering kinetics Therefore, Eq. (|21jl can be viewed 
as a generalization of the Cahn-Hilliard equation for long-range interactions j^ . 

We shall now assume that the potential V(4>) is a symmetric function of <p with two 
minima at ±u such that V(±u) = 0. This is a typical prediction of the statistical theory 
in the two-levels approximation where C(<p) = (1/2){(1 + <fi) ln(l + <fi) + (1 — <p) ln(l — <p)} 
after normalization 30] . However, our goal here is to construct more general solutions of the 
Q.G. equations enjoying the same properties. The potential V{<f>) will have two minima if 
the equation C'(u) = —R 2 f3u has non zero solutions. For typical situations, this will be the 
case if j3 > f3 c = —C"(0)/R 2 . The existence of a critical (negative) temperature is of course 
reminiscent of a second order phase transition in thermodynamics. Note that the functional 
(I26j) considered by Turkington and collaborators does not satisfy this condition so that the 
limit of small Rossby radius R — > is different in that case and does not lead to a thin jet 
structure. 

From now on, we assume that j3 > (3 C . The stationary solutions of the Cahn-Hilliard 
equation (}3T]) satisfy the differential equation 

R 2 A(j) = V\<j>) - h + a, (32) 

where a is an integration constant. The same equation is obtained by minimizing the Landau 
free energy ()30|) at fixed V. For R — > 0, these solutions describe "domain walls" in the theory 
of phase ordering kinetics. They connect regions of uniform potential vorticity q = B ±u 
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separated by a wall. This is precisely the structure of Jupiter's great red spot where the 
PV gradient is concentrated in a thin annular jet scaling with the Rossby radius. This 



analogy has been developed in 



28, 



30] and it was shown to give a good agreement with the 
GRS. The present approach makes a connexion with the theory of phase ordering kinetics 
and Cahn-Hilliard equations and is valid for an arbitrary jy-shape potential V(q). On a 
technical point of view, we differ from [s(3] by treating the free energy as a functional of q 
(instead of ?/>), using the short-range expansion (}2T?|) of ip. This makes a close link between 
the free energy ()28|) and the Landau free energy f)30j) . Having established this link, we can 
now directly use standard results of domain wall theory [8| to describe the great red spot. 
Treating the topography as a perturbation of order R, we find that the wall profile is the 
solution of the equation 



2, 



r 2 ^ = n*\ (33) 

with boundary conditions 0(±oo) = ±u, where £ is a coordinate normal to the wall. This 
is like the equation of motion for a fiducial particle in a potential V. The first integral is 



Rd<p/d£ = y/2V((f>). This result can be used to give the surface tension 

o = R 2 J + {^\ d^ = Rj + v/2V(0)#. (34) 

Finally, linearizing Eq. (|33|) around <fi = ±u gives u =F ~ exp [-^V"{u)i/R] so that 
the typical extension of the PV wall is L = Rj \/V"(u). The domain wall curvature due 
to the topography h{y) can be obtained as follows. Close to the interface, V0 = difi/d^n 
where n is a unit vector normal to the wall. Introducing the curvature K = V ■ n, we get 
A(f) = Kd(j)/d£, + d 2 4>/d^ 2 . Therefore, Eq. becomes to first order 

R2 ( K ^ + ^) =n*) (35) 

Multiplying by dip/d^ and integrating over the wall, we obtain the "curvature-topography" 
relation 

—?L = ai -h(y). (36) 

This relation determines the elongation of the GRS under the effect of an underlying topog- 
raphy (more precisely a deep shear layer). The above equations are similar to those derived 
by Bouchet & Sommeria [3i| except that they apply to the PV wall instead of the velocity 
jet. They are also expressed in the case of an arbitrary potential V(<p) so that they can be 
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used to construct a wider class of models. By construction, all these models are granted 
to be nonlinearly dynamically stable via the Q.G. equations. It has already been shown by 
Bouchet & Sommeria |3(| that this description gives a fair agreement with the structure of 
jovian vortices when S[q] is the Fermi-Dirac entropy, corresponding to the two- levels ap- 
proximation of the statistical theory. However, some discrepencies have also been noted (in 
particular, the predicted surrounded shear is 3 times smaller than its real value). It would be 
interesting to determine whether an equally good (or even better) agreement can be achieved 
by other stationary solutions of the Q.G. equations, that do not necessarily correspond to 
statistical equilibria. This could be used to test the power of prediction of the statistical 
theory and determine a class of relevant functions C(q) for jovian flows. The present paper 
has given the theoretical tools for constructing numerically such general solutions and this 
problem will be considered elsewhere. These methods can also be extended to the more 
realistic shallow- water (SW) equations 2^|. 

V. CONCLUSION 

In this paper, we have proposed a parametrization of 2D geophysical flows in the form 
of a generalized Fokker-Planck equation and we have showed the connexion with Cahn- 
Hilliard equations in the limit of small Rossby deformation radius, leading to short-ranged 
interactions. This equation is associated with a generalized entropy functional which is fixed 
by a prior vorticity distribution in the sense of Ellis et al. | . This equation can thus describe 
the evolution of large-scale motion when the small-scale forcing has established a permanent 
vorticity distribution. Its domain of validity is therefore limited close to equilibrium like in 
Onsager's linear thermodynamics. However, it may remain of practical interest even if we 
are far from equilibrium. The parametrization depends on a prior vorticity distribution (or 
on a generalized entropy) which is unknown in general and which must be adapted to the 
situation. Several prior distributions/entropies that fall in the same "class of equivalence" 
l7| should give similar results. We have shown, however, that the Fermi-Dirac entropy used 

n n 

by Bouchet & Sommeria [30( and the entropy used by Ellis et al. |6[ give different results in 
the limit of small Rossby radius. They thus belong to different "classes" . The Fermi-Dirac 
entropy may be more appropriate for describing the Great Red Spot but we have proposed 
that it may correspond to a prior vorticity distribution created by a small-scale forcing with 
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two intense peaks rather than to a free evolution of the system with two vorticity levels 
as in jscj. These results can be extended to the more realistic shallow- water equations so 
their domain of application is wide. The challenge now amounts to finding relevant forms 
of prior vorticity distribution for each specific situation. This could be achieved by writing 
stochastic processes for the generation of potential vorticity (in preparation). 

Of course, this thermodynamical parametrization assumes that the statistical theory 
works well. However, it has been realized in many occasions that the ergodic hypothesis 
which sustains the statistical theory is not fulfilled everywhere so that the statistical predic- 
tion (assuming efficient mixing) is not truly reliable. Indeed, the system can be t rapped in a 
stationary solution of the 2D Euler equation which is not the most mixed state P. fiol. I20L I33I] . 
In general, 2D vortices (and galaxies in astrophysics) are more confined than predicted by 
the statistical mechanics of violent relaxation [13]. The system tends toward the statistical 
equilibrium state during violent relaxation but cannot attain it: the fluctuations die away 
before the relaxation process is complete. This phenomenon is referred to as incomplete 
relaxation [l]]. It is an important obstacle for the general application of the statistical 
mechanics of violent relaxation. 

This effect of incomplete relaxation can be taken into account in the relaxation equations 
by using a space dependent diffusion coefficient of the form (|2*1|) which is related to the 



local fluctuations of the vorticity 







24j . This can freeze the system is a sub-domain 



of space, in a sort of "maximum entropy bubble" |3J| surrounded by an un-mixed region 
q ~ which is poorly sampled by the flow. This justifies dynamically why the statistical 
equilibrium state is not always reached in practice. This is an interesting property of the 
relaxation equations which cannot be obtained with the Turkington-Whitaker algorithm 
27| as it assumes complete relaxation. Alternatively, we can try to construct nonlinearly 
dynamically stable stationary solutions of the Euler equations that correspond to confined 
structures. In that case, the relaxation equation (|21j) can be used as a numerical algorithm 
to construct solutions that are not catched by the statistical theory. In that context, the 
functional S[q] is interpreted as a H-function. For example, Tsallis functional is a particular 
H-function leading to "polytropic vortices" with a compact support (the enstrophy is a 
particular case corresponding to q = 2) [7]. Their confinement is interpreted as an effect 
of incomplete relaxation due to lack of mixing/ ergodicity. In case of incomplete relaxation, 
the goal is not to predict the metaequilibrium state (this is probably an impossible task) 
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but simply to devise a general method for constructing robust stationary solutions of the 
2D Euler equations in an attempt to reproduce observed phenomena. This is already a 
non-trivial problem by itself. Of course, the statistical theory can be a valuable guide to 
select physically motivated solutions but the relaxation equation (|21j) can be used to depart 
from the statistical prediction and construct a larger class of models. The same ideas can 
)e developed for collisionless stellar systems and other systems with long-range interactions 

The relaxation equation (|2ip providing a parametrization of geophysical flows has been 
justified from a method similar to Onsager's linear thermodynamics or from a Maximum 
Entropy Production Principle [2 , 23 1 . It would be nice to have a more rigorous justification 
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22| using 



from "first principles" by using kinetic theories. An attempt has been made in 
formal methods. However, the domain of application of these methods is limited and does 
not correspond to the very nonlinear stages of "violent relaxation" that we would like to 
describe in concrete situations. This is why more phenomenological approaches are prefered 
until we have a more satisfying theory. Still, the kinetic theory tends to confirm the general 
structure (drift-diffusion) of the relaxation equations obtained with the MEPP. Finally, we 
note that generalized Fokker-Planck equations of the form (j21|) appear in other domains 
of physics (porous media, Brownian particles in interaction,...) and biology (chemotactic 
aggregation of bacterial populations) where they have a different status and a different 
interpretation 0, Therefore, the general study of these equations is of considerable 
interest, independently of the context of 2D hydrodynamics. A systematic study of this 
class of equations has been undertaken in 



APPENDIX A: RELATIVE ENTROPY 



In this Appendix, we show that the relative entropy ()10|) can be seen as a Legendre 
transform of the mixing entropy (jSJ) when the constraints on r n> i are treated canonically. 
The variational principle leading to the Gibbs state © can be written 

5S - a JT n - (35E -a5T- I Q{v)8 ( j pdaj d 2 r = 0, (Al) 

n>l \J / 

where we have distinguished the robust constraints E, T and the fragile constraints r n >i- 
In the point of view of , the moments r n >i are treated microcanonically and we must 
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ultimately relate the Lagrange multipliers a n to the constraints. In the point of view of 
0| , these constraints are treated canonically by fixing the Lagrange multipliers a n . If we 
regard the PV levels as different species of particles, this is equivalent to fixing the chemical 
potentials instead of the total number of particles in each species. We are led therefore to 
define a relative entropy 

S x = S-J2 a nT n . (A2) 

n>l 

This is similar to the Legendre transform F = E — TS in usual thermodynamics when we 
pass from the entropy (in a microcanonical description where the energy is fixed) to the free 
energy (in a canonical description where the temperature is fixed). Explicitly, 



S x = - p hip + a n° n 
J L n>l 

Introducing the prior vorticity distribution 



d rdcr. 



(A3) 



(A4) 



we get 



S x = - \ P ln 



d 2 rdcr, 



(A5) 



which coincides with Eq. (1TU|) if we set /i = p/x{&)- 



APPENDIX B: DIFFUSION COEFFICIENT IN THE Q.G. MODEL 



In this Appendix, we derive an expression for the diffusion coefficient appearing in the 



parametrization 
Robert & Rosier 



Q) related to the Q.G. model. We extend the arguments developed by 
24| and Chavanis et al. [12] 



(see also 

diffusion coefficient can be estimated by the formula 

1 



2l| ) for the 2D Euler equation. The 



D 



4 



ru 2 (r, t) 



(Bl) 



where r is the decorrelation time of the system. Equation (jBl|) corresponds to the general 
Taylor expression of the turbulent viscosity. In the Q.G. model, the velocity is related to 
the potential vorticity by 



u(r,t) = J q(r',t)K{r-r')d 2 r' 



(B2) 
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K(r - r') = z x J-K x (^-^) t^— (B3) 



with 

K(r - r') = 7, x 

27ri? V i? y [r — r' 
where ifi(x) = —K' (x) is a modified Bessel function. Therefore, the fluctuations of velocity 
are induced by the fluctuations of potential vorticity so that the expression of the diffusion 
coefficient becomes 



D = ~r J K(r - r') ■ K(r - r")g(r', t)q(r", t)d 2 r'd 2 r" 



(B4) 



We shall now consider that the scale of the spatial correlations e is smaller than the Rossby 
radius R and write the correlation function as 



g(r', t)q(r", t) = e 2 g 2 (r', t)5(r' - r"). (B5) 
Substituting this expression in Eq. (IB4I) . we obtain 



D = TT 1 ^ / ^if^^ )^',t)d 2 v'. (BO) 



16^2 J \ R 

Making a local approximation, we can rewrite the foregoing expression in the form 



re 2 r +oc 

D = —q 2 (r,t) JC?(£)£de, (B7) 
8vr J e/R 



where q 2 = q 2 = (q — q) 2 is the local centered PV enstrophy. The integral can be approxi- 
mated by ln(i?/e) so we obtain 

D = ^\n(-)q 2 (r,t). (B8) 



8tt V e , 

We see that the introduction of a finite Rossby radius regularizes the logarithmic divergence 
at large scales encountered in the case of the 2D Euler equation [12J, 
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24]. In the present 

case, the integral has to be cut-off at R, the typical range of the interactions, while in the 
former case it had to be cut-off at the typical vortex or system size. Finally, the decorrelation 
time can be estimated by r ~ e 2 /D. This yields 

D = Ke 2 ql'\ (B9) 

where K = h\{R/e)} 1 ^ 2 is a constant of order unity. 
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APPENDIX C: NON-EQUILIBRIUM DISTRIBUTIONS 



The relaxation equation (|21|) gives the evolution of the coarse-grained PV q(r, t) out of 
equilibrium. This evolution respects the conservation of energy and circulation and depends 
on the prior x( a ) through the function C(q). We may wonder about the expression of the de- 
tailed PV distribution p(r, a, t) out of equilibrium. One idea is to determine this distribution 
by using a maximum entropy principle. Specifically, one may argue that p(r, a, t) maximizes 
the relative entropy ()10j) at fixed coarse-grained PV q(r, t) = J pada and normalization 
J pda = 1. This yields 

p(r,v,t) = ^j )X (a)e-^\ (CI) 
where Z and $ are Lagrange multipliers determined by 

Z(r,t) = J X (a)e-^da, (C2) 

?(r, t) = -^—^ I X {a)ae-^da. (C3) 

This is similar to the equilibrium problem with $(r, t) instead of (3ip + a. In particular, the 
relation (|C3J) can be rewritten 

g(r,t) = F[$(r,t)], (C4) 

where F is the same function as in Eq. (JJJ), entirely determined by the prior x{ a )- Thus, 
knowing q(r, t) from the evolution equation ([21)1. we can inverse (at each time) Eq. (|C4|) to 
obtain $(r, t), then p(r,a,t). On the other hand, from the relations q'[§(r,t)} = —q2{r,t) 
and C'[q(r,t)] = -$(r,t) leading to q'[®(r,t)) = -1/C"[q(r,t)), we get 

'^• t) = a4^yy (C5) 



Therefore, this approach justifies, or at least is consistent, with the use of Eq. (|23|) out of 
equilibrium. 
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